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In this paper, we view fluctuating fronts made of particles on a one- dimensional lattice as an 
extreme value problem. The idea is to denote the conflguration for a single front realization at time 
t by the set of co-ordinates {kiit)} = [fci (t), fc2(t), . . . , fcjv(t) (*)] of the constituent particles, where 
N{t) is the total number of particles in that realization at time t. When {ki{t)} are arranged in 
the ascending order of magnitudes, the instantaneous front position can be denoted by the location 
of the rightmost particle, i.e., by the extremal value kf{t) = max[fci (i), feCt), . . . , fcjv(t) (i)]- Due to 
interparticle interactions, {ki{t)} at two different times for a single front realization are naturally 
not independent of each other, and thus the probability distribution Pk^ (t) [based on an ensemble 
of such front realizations] describes extreme value statistics for a set of correlated random variables. 
In view of the fact that exact results for correlated extreme value statistics are rather rare, here 
we show that for a fermionic front model in a reaction-diffusion system, Pfc^. (t) is Gaussian. In a 
bosonic front model however, we observe small deviations from the Gaussian. 

PACS numbers: 02.50.-r, 02.50.Ey, 45.70.Qj 



I. INTRODUCTION 

Extreme value statistics of random variablesplays a 
diverse role in physics, chemistry and biology 0, S H- 
The topic concerns the probabihty distributions of the ex- 
trema (i.e., the maximum kmax or the minimum kmin) of 
a set of N random variables {ki, k2, ■ ■ ■ , k^} in the hmit 
N oo. When the random variables ki are uncorre- 
lated, the probability distribution of kmin belongs to one 
of the three universality classes jj| , but the identification 
of similar universality classes for the extreme value statis- 
tics of correlated random variables is still largely an open 
problem. A few results relating to extreme value statis- 
tics for correlated random variables in physics, computer 
science and mathematics have been obtained in the re- 
cent past nevertheless, any exact result that can be 
obtained for correlated random variables is an important 
addition to the present state of knowledge. 

From this perspective, in this paper, we present two 
main results relating to fluctuating fronts made of dis- 
crete particles on a one-dimensional lattice. Before we 
proceed further with our formulation of the problem, 
we must note that an intriguing connection between the 
extreme value statistics of correlated random variables 
and travelling fronts have already emerged from the re- 
cent works 00- To be more precise, these works have 
demonstrated, for the models they studied, that the cu- 
mulative probability distributions of extreme values for 
correlated random variables admit propagating front so- 
lutions, wherein the variance of the extremal variable 
is the front width itself. Our formulation here, how- 
ever, is completely the other way round: namely that 
our systems consist of many interacting particles, where 
the dynamics of the systems already admits front solu- 
tions propagating into unstable states. Although in a 
deterministic mean-field description, these fronts propa- 



gate with a fixed speed and a fixed shape at long times, 
due to the presence of stochasticity involving many par- 
ticles, the front in a given realization of the system does 
not move with a uniform speed even at long times — 
instead, the front speed averaged over an ensemble of 
front realizations approaches a constant in time at long 
times. Moreover, as a result of the inherent stochasticity 
in these systems, the individual front realizations that 
are initially aligned with each other do not remain so at 
a later time; instead their displacement w.r.t. each other 
keeps increasing with time (see Fig. 4 of Ref. Q for an 
illustration). As explained below, it is the dynamics of 
the individual front realizations in the ensemble that we 
pose as a correlated extreme value problem in this paper. 

The correspondence between the extremal value statis- 
tics and the fluctuating fronts in these systems is easily 
made by first noticing that in any realization of these 
systems, the front position can be denoted by the instan- 
taneous position kf{t) of the foremost (or the rightmost) 
particle li|l3- The interest then lies in the probability 
distribution Pk^ (t) , which describes the statistics of the 
front position in time for an ensemble of front realiza- 
tions. Secondly, in a snapshot of one single realization, 
the configuration of the system is described by the lo- 
cations of the particles (as random variables) {ki{t)} = 
[ki{t), k2{t), . . . , A:jv(t)(i)], where N{t) is the total number 
of particles in that realization at time t. Then the instan- 
taneous front position kf(t) in this formulation is then 
simply the extremal value max[fci(t), k2{t), . . . , fc7v(t)(i)]- 
Due to the interparticle interaction within the system 
defined by the microscopic rules of the dynamics, {ki{t)} 
are naturally not independent of each other, and thus 
Pkf (t) simply describes the statistics of the extreme for 
a set of correlated random variables. 

In this paper, we consider two different systems that 
admit front solutions propagating into unstable states: 
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fa) a fermionic reaction-diffusion system A^A+A 0, 
in Sec. IHI where we show that Pk^ (t) is Gaussian, 
and (b) the so-called (bosonic) clock model in Sec. 
mil where Pk^ (t) has small deviations from the Gaus- 
sian. It is important to note here that the front solutions 
in these models have been analyzed before, in the sense 
that both the front speed v = limt^oo d{kf{t)) /dt and 
the front diffusion coefficient Df = limt^ao d{[kf{t) — 
vtY)/dt, respectively based on the first and the second 
moments of Pkf(t), have previously been analyzed and 
numerically measured 0, [li IH El [H [l| . The higher 
(than second) moment of Pkf{t), or Pk;{t) itself, how- 
ever, have not been extensively studied before. 

The paper is finally ended with a discussion in Sec. IIVI 



II. A FERMIONIC REACTION-DIFFUSION 
MODEL AND GAUSSIAN BEHAVIOUR OF Pfe^(t) 

In this model, we consider a one-dimensional lattice on 
which at most one A particle is allowed per lattice site at 
any instant — hence the model is named fermionic. The 
particles can undergo the following three basic moves, 
shown in Fig. ^ (i) A particle can diffuse to any one of its 
neighbour lattice sites with a diffusion rate -D, provided 
this neighbouring site is empty, (ii) Any particle can give 
birth to another one on any one of its empty neighbour 
lattice site with a birth rate e. (iii) Any one of two A 
particles belonging to two neighbouring filled lattice sites 
can get annihilated with a death rate W . 
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FIG. 1: The microscopic processes tliat take place inside the 
system: (i) A diffusive hop with rate D to a neighboring 
empty site; (ii) Creation of a new particle on a site neigh- 
boring an occupied site with rate e; (iii) Annihilation of a 
particle on a site adjacent to an occupied site at a rate W . 

The lattice indexed by k that we consider in this prob- 
lem is semi-infinite. The left boundary is impenetrable 
— no particle can diffuse across the left boundary located 
on the left of the lattice site A: = 0, while the system is 
of infinite extent on the right side. Following the usual 
convention, we start with a step initial condition, i.e., 
at time f = 0, there exists a finite k-^igtAi such that all 
lattice sites < fc < fcright are occupied and k > bright 



are empty. This system then admits a fluctuating (and 
propagating) front solution for t > 0. 

Earlier work on models of this type has appeared in 
Refs. mill 113,155. In the general" case there are essen- 
tially only two nontrivial parameters in our model, e.g., 
the ratios D/e and D/W, since an overall multiplicative 
factor simply sets the time scale. When these ratios tend 
to infinity, the front speed approaches its mean field value 

m- 

For an ensemble of front realizations, let us denote the 
probability distribution for the foremost occupied lattice 
site to be at lattice site kf at time t by Pkj,{t). The 
evolution of Pk^ (t) is then described by 



dt 



£) Pkf-l 



[D + e) Pk, 



r) pcmpty , rrr pocc 

^pompty ^ pocc 



(1) 



Here -Pfc™(t) and P^J^^^^it) respectively denote the joint 
probabilities that the foremost particle is at site kf and 
that the site kf — 1 is occupied and empty. Clearly, 

Pkfit) = P^rii) + PkT'(i)^ and Efc, PkAt) = 1- The 
first term on the r.h.s. of Eq. describes the increase in 
Pkf (t) due to the advance of a foremost occupied lattice 
site from position kf — 1, while the second term describes 
the increase in Pkf {t) due to the retreat of a foremost oc- 
cupied lattice site from position kf+l. The third and the 
fourth terms respectively describe the decrease in Pkf (t) 
due to the advance and retreat of a foremost occupied 
lattice site from position kf. It is clear from this formu- 
lation that the dynamics of Pkf (i) is effectively obtained 
only from the coupled interaction between the foremost 
particle and the site just behind it. 
In addition to Eq. we have 



P^;%t) = pkfMt)Pkfit), 



(2) 



where pkf-i{t) is the conditional probability of having 
the (A;/ — l)th lattice site occupied. At large t, pkf-i{t) 
should be independent of kf and t, and one can replace 
Pkf-i(t) by p in Eq. (0), where the numerical value of 
p depends only on those of Z), e and W . Similarly, the 
set of (time and ^/-independent) conditional occupation 
densities Pkf-mit) for m > 1 can be thought of as de- 
termining the front profile in a frame moving with the 
foremost particle of each front realization (see Fig. 5 of 
Ref. 10] for an illustration). 

With the condition Pkf it) = Pk'^'it) + Pk"''^^^{t), and 
the notation q — p{W ~ D), at large t, Eq. |^ can be 
rewritten as 

^^l{2D + e + q) [Pkf+, + Pkf-r - 2P.,] 

-\[e~q) [Pkf+i-Pkf-i] , (3) 

which is clearly a diffusion equation for Pkf (t) with a 
drift. After having aligned the locations of the fore- 
most particles for all realizations in the ensemble, say 
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at kf = kin at time Un > 1 [i.e., Pk^it-m) = 6kf,kiJ' "^^ 
are interested in the solution of Pkf{t)- In fact, Eq. Q 
can be solved by taking a discrete Fourier transform in 
fc/, but due to the redundancy of the wavevector modulo 
any multiple of 27r, the magnitude of the wavevector has 
to be kept confined only within the first Brillouin zone 
(— 7r,7r]. Then for At = t — ti^ ^ 1, it is easily seen 
that the dominant contribution to Pkf (t) comes from the 
wavevector in the neighbourhood of zero, yielding |16| 



exp 



{kf - fcin - vAt)"' 

iDfAt 



(4) 



Here, v = e ~ q is the front speed and -D/ = "20 + e + q 
is the front diffusion coefficient, as already derived as the 
first and the second moment of Pk^ (t) in Ref. ■ 



III. THE CLOCK MODEL AND THE 
NON-GAUSSIAN BEHAVIOUR OF Pk^ (t) 

The clock model was originally invented in the con- 
text of the largest Lyapunov exponent for a gas of hard 
spheres this model, one considers a system of 

N clocks with integer readings {fci}. The dynamics of 
the clocks involve binary "collisions" between any two 
randomly chosen clocks in continuous time. In a colli- 
sion between two clocks with pre-collisional readings ki 
and kj, the post-coUisional readings of both clocks are 
updated to max[/ci, kj] + 1. 

In the clock reading space, which can be imagined 
as a one-dimensional lattice, the number of clocks Nk 
with readings k or higher for any realization of the clock 
model admits a fluctuating (and propagating) front solu- 
tion Clock model allows more than one clock with 
the same reading and hence the model is bosonic. Con- 
ventionally, all clock readings in any realization are ini- 
tially (i.e., at t = 0), taken to be equal to zero — for 
the propagating front, this corresponds to the step ini- 
tial condition 1131. Once again we denote the largest 
clock reading in any realization at time t by kf{t). 

In the deterministic mean- field limit, the prop agating 
front in the clock model is a pulled front and if 

the time is rescaled in order to have the mean collision 
frequency of a single clock equal to unity, the front prop- 
agates with a speed v* = 4.31107 . . . (13]. However, due 
to stochasticity effects associated with discreteness effects 
of the clocks and their readings, in the limit of asymp- 
totically large TV, the front speed v and front diffusion 
coefficient Df, which could be measured following the 
procedure described in the last paragraph of Sec. ^ have 
the property that (v* — v)oc l/ln^TV and £»/ oc 1/ln^iV 
[T^ . Thus, the clock model is an example of a fluctuating 
"pulled" front HHIll. 

To write a master equation for Pkj, {t) defined over an 
ensemble in the clock model, it may be argued that the 
reading of any clock in any realization can only increase 



with time; and thus, Pk^ (t) can increase when in a real- 
ization, one of the clocks with largest reading kf — 1 is 
involved in a collision with another one. Similarly, Pk^ (t) 
can decrease when in a realization, one of the clocks with 
largest reading kf is involved in a collision with another 
one. If we now denote the conditional probability of the 
number of clocks with largest reading k / to be nkj (t) at 
time t by V{nkf , t), the master equation for Pkj (t) reads 



dPk, 
dt 



C{7ikf-i)V{nkf-i,t) 



Pk 



Y.C{nk,)V{nk,,t) 



(5) 



Here, C{nkj) is the rate of collisions that involve a 
clock with reading kf for a realization with kf as the 
largest of the clock readings. From Eq. one 
might now further argue that at large the quanti- 
ties within the large square brackets in Eq. (|SJl are 
independent of t and fc/, and thus at large i, Eq. (O 

dPk, _, 

should reduce to a form — 7— = c[Pkf-i — Pkfli where 



dt 



J2C{nkf)V{nkf,t) 



at large t. However, for 



any finite value of N, the simple-minded replacement of 



J2C{nkf)V{nk,,t) 



by a t and fcf -independent quan- 



tity c in Eq. lO at large t is incorrect for the clock model 
— caused by the fact that Vi rik, , t) does not become 
time-independent at large t |20j | — as we argue below. 

The observation we make, in order to argue that 
V{nkj,t) does not become independent of t at large t, 
is that the largest clock reading in any realization does 
not increase smoothly in time with a rate v even at large 
t. Instead, after attaining a new integer value, the largest 
of the clock readings for any given realization does not 
change for some time-interval (hereafter denoted by St) 
of typical magnitude 1/v before attaining the next inte- 
ger value 10. Generally speaking, during any of these 
time- intervals, the number of clocks with the largest clock 
reading in any realization increases with time; and the 
number of clocks with the largest reading at any instant 
in a given realization depends on how long the largest 
clock reading remains unchanged at its value. The con- 
ditional probability P{nkf, t) can thus be written as 



P{nkf,t) 



d{St) pi{St,t) p2{nkf,St,t) 



(6) 



With Eq. © in the back of our minds, we now return 
to the statement to the second sentence of Ref. [Tsj : 
namely that front propagation in any realization of the 
clock model is coded in the sequential values of the time- 
intervals {5ti} between the consecutive changes of the 
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largest clock reading. In this description, the point to 
note is that the St^ values are very strongly correlated 
with each other; e.g., a large St is almost always fol- 
lowed by several small values of 6t and vice versa (the 
large or smallness of 6t are decided in comparison to 1/v) 
d, 13 ■ Due to such strong dependence of of the St values 
on the evolution histories of individual realizations, it is 
easily conceivable that the shape of the probability dis- 
tribution pi{St,t) lacks ^-independence at large t. where 



1710 




i ' ' ' ' ' ' ' ' 

100 200 300 400 

FIG. 2: Top figure: Simulation results for {kf{t) — kin) as 
a function of t — tin- Bottom figure: Simulation results for 
{[kf(t) — kin]'^) as a function of t ~ tin- Apart from an initial 
transient for t — tin < 10, {[kf{t) — kin]'^) increases linearly with 
t, indicating that the front wandering is diff'usive. Crosses 
correspond to N = 10"* (front speed v = 4.08 and front dif- 
fusion coefficient Df = 0.112) and filled circles correspond to 
A'^ — 10'' (front speed v — 4.17 and front diffusion coefficient 
Df = 0.056) in both figures. 

pi{St,t) is the probability that the largest clock reading 
became fc/ at time (t — 6t) and remains so until time t, 
and p2{nkj , St) is the probability of having rikf clocks at 
time t when the largest clock reading became kf at time 
(t — 6t) and remains so until time t. Using Eq. lO, the 
t-dependence of V{nkj , t) can then be argued in terms of 
the t-dependences of pi{5t,t) and p2{nkf ,St,t). 



The t-dependence of p2{nkf,St,t) can be argued in a 
similar way. In realizations for which the largest clock 
value became fcj at time (t — St) and remains so until 
time t, how many clocks share the largest clock reading 
at time t depends on the time-dependence of the num- 
ber of clocks nkj-i with clock readings (fc/ — 1) between 
times (t — 6t) and t — after all, any clock that attains a 
reading kf must come out of a collision that involves a 
clock with reading {kf — 1). Between times {t — St) and 
t, — 1 changes also with time, and thus the probabil- 
ity distribution p2{nkf,St,t) inherently connects to the 

fuctuations in the shapes of individual front realizations 
LH^. These fluctuations have a typical correlation time 
oc In^ dlllll^. For A^ ^ oo, this correlation time 
also becomes large, and one therefore expects the shape 
of p2{nkj , St, t) to also depend on t via the strong depen- 
dence of Ukf-i on the evolution histories of individual 
realizations at earlier times. 



0.45 
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FIG. 3: Pkf (t) for the clock model; pluses: data for A^ = 10* 
at t = tin + 99, circles: data for A^ = lO" at t = tin + 297, 
crosses: data for TV = 10^ at t — tin + 198 and diamonds: data 
for N = 10® at t — tin + 396. Solid line: normalized Gaussian 
distribution with mean zero and variance unity. Dotted line: 
numerically obtained curve for the collapsed data. 

With no further simplification of Eq. Q possible, let 
alone an exact solution for Pk^ {t) like in Eq. Q), we can 
only study Pk^ (t) for the clock model only via simulation. 
Our (molecular dynamics) simulation methods are as fol- 
lows: we choose an ensemble of A/" = 50000 realizations 
of A^ = 10* clocks and set all clock values zero at t = 0. 
We then let each realization evolve until time tin — 800 
units. At tin, we align the different realizations in such 
a way that the largest of all the clock values coincide at 
k = kin- We then follow the locations of the largest clock 
values for each realization until t — tin — 400. We also 
repeat the calculations for the same values of A/", tin and 
kin but for A^ = 10^. The ensemble average {kf{t) — kin) 



5 



and {[kf{t)-kin]'^) for A/" = 50000 as a function of (t-tin) 
both for TV = 10^ and 10^ have been shown in Fig. |21 

To obtain Pk^ (t) numerically from the above data, we 
now proceed in the following way. First we select two 
different time instants for each value of N to take snap- 
shots of the entire ensemble of fc/-values: for N — 10^, 
we choose t = + 99 and tin + 297, and for = lO*^, 
we choose t = tin + 198 and tin + 396. Having used 
the best fit method from the data of Fig. [S] we then 
identify the location of the mean front position fco(t) 
and the standard deviation CTo(i) [coC^) effectively be- 
haves as ~ y/2Df{t — tin) as seen in the bottom plot 
of Fig. 121 for large (t — tin)], for two different values of 
N at these different time instants. Finally, with the 
histograms /^{kf ,t) /[Afao{t)] plotted as a function of 
[kf — ko{t)]/ao{t), where M{kf,t) is the number of re- 
alizations with largest of the clock values kf aX time t, 
we expect a good data collapse, and the corresponding 
curve then gives us the normalized Pkj (t) . Notice that 
the procedure that we followed to obtain A:o(t) and cro(t) 
[and subsequently the numerical curve for Pkj (t)] at the 
above time instants does not guarantee {kf{t))—ko{t) = 
and ([/i:/(t)-fco(t)]2) = crg(t); instead, the (A;/(t)) - fco(t) 
and the {[kf{t) — k^ity^) j a^if) values are in fact very 
close to zero and unity respectively. 

This data collapse is shown by means of the numeri- 
cally obtained dotted curve in Fig. |31 Further analysis of 
the data ( not presented here ) clearly shows that the dot- 
ted curve does not belong to any of the known universal- 
ity classes for the extreme value statistics of uncorrelated 
random variables 0/; instead, it appears to resemble the 
normalized Gaussian distribution rather closely. To fa- 
cilitate comparison, we therefore plot Pj.^ (t) against the 
normalized Gaussian distribution (with mean zero and 
variance unity). It is clear from Fig. |3| that the dot- 
ted curve is positively skewed; direct measurement of the 
third cumulant from the data also confirms this positive 
skewness behaviour of Pkf (t)- The most noteworthy fea- 
ture is the longer right tail of the collapsed data than the 
left tail, implying that the probability for large positive 
deviation around the mean for the clock model is larger 
than that of large negative deviation. This is indeed con- 
sistent with positively skewed Pkf (t) — as stated before, 
(fc/(t) — ko{t)) ~ for all snapshots. 

While Fig. O certainly provides an example of devia- 
tion from Gaussian statistics when the fluctuating front 
propagation is seen as a correlated extreme value prob- 
lem, it also provides an interesting perspective from the 
point of view of fluctuating front propagation literature. 
As already mentioned before, clock model is an exam- 
ple of fluctuating "pulled" fronts, and the expression for 
V and the scaling for Df due to the discrete particle 
stochasticity effects in the limit of asymptotically large 
values of N are known for the last few years. It is also 
known that over a time interval At at large t, the second 
moment of Pfe^(t), i.e., ([fc/(At) - fcin - wAt]^) ^2DfM 
for all values of N . Figure |3| however shows that the 
information regarding the second moment is clearly not 



enough to characterize Pkj {t). Nevertheless, the data 
collapse shows that at large t, Pkj{t) = P[{kf — fcin — 
vAt)/ y/2DjAi]/ y/2DjAi (i.e., the dotted line in Fig. 
O is a characteristic curve for the clock model, and this 
characteristic curve is not Gaussian for the values of N 
studied here. The statement that "the front wandering 
is diffusive" at any value of N for the clock model, there- 
fore, has to be interpreted only in the sense that the 
second moment of Pkf (t) increases linearly with time at 
large t for any value of N. 

Whether the deviation of the dotted line from the 
Gaussian is due to the fact that we have not used ex- 
tremely large values of N is however not clear. It is well 
known that to observe the 1/ln^iV scaling of {v* — v) 
and the l/\n'^ N scaling of Df for fluctuating "pulled" 
fronts one needs to take N extremely high [ij, [lj| . Di- 
rect molecular dynamics simulations of the clock model 
for N > 10^ are prohibitively slow. The existing simu- 
lation methods at much higher values of TV are not only 
quite intricate, but they also do not follow the exact dy- 
namics of the model for all clocks. This particular point, 
therefore, is left here for further investigation in future. 



IV. DISCUSSION 

In this paper, we have analyzed front propagation in 
discrete particle systems on a one-dimensional lattice as 
extreme value problems. In these systems, the positions 
of the particles can be thought of as random variables, 
and these random variables under consideration are obvi- 
ously strongly correlated with each other. We have seen 
that in the case of the fermionic reaction-diffusion model, 
the extreme value problem follows Gaussian statistics. It 
clearly does not belong to any of the classes pertaining to 
extreme value statistics of uncorrelated random variables. 
For the (bosonic) clock model however, we see that the 
extreme value statistics has a small deviation from the 
Gaussian, and additional analysis (now presented here) 
also clearly shows that the probability distribution does 
not belong to any of the known universality classes for the 
extreme valus statistics of uncorrelated random variables. 
However, due to the unavailability of any analytical tool, 
the characterization of this distribution has proved elu- 
sive. Whether the small deviation from the Gaussian is 
caused by the fact that we have not used extremely high 
values of iV for our simulations thus remains an open 
question. 
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